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Relating B s Mixing and B s — > with New Physics - An Update 

E. Golowich 
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This document describes my talk (based on work by JoAnne Hewett, Sandip Pakvasa, Alexey Petrov, Gagik 
Yeghiyan and myself) given at the 2011 Meeting of the Division of Particles and Fields of the American Physical 
Society (8/9/11-8/13/11) hosted by the Physics Department at Brown University. We perform a study of the 
Standard Model (SM) fit to the mixing quantity AMb s in order to bound contributions of New Physics to B s 



(N 
> 



(NP). In some cases, this constrains NP amplitudes for B s — > M + M to lie below the SM component. 



Oh 

1. Introduction 

Here, I describe a calculation [1] carried out by JoAnne Hewett, Sandip Pakvasa, Alexey Petrov, Gagik 
Yeghiyan and myself (hereafter denoted collectively as GHPPY) in which we study the transition B s — > /J + /i~ 
O ■ in the context of both the Standard Model (SM) and New Physics (NP). As described in Section 2, our approach 
involves the mixing mass difference AMb s as well. The 'Update' mentioned in the title refers to two items. 
The first is to apply to our analysis the LHC bound for the branching ratio BrB g ^p 1 + ^- which was recently 
made public by Ref. [2]. For example, this affects the exclusion region in parameter spaces of NP models (cf 
i— i, Figs. 3,5). The second update refers to a change in one of the inputs in the published article of Ref. [1] which 
modifies our determination of AAIb 3 (c/Eq. (4)) and some other some of the numerical results displayed in 
this report. This is also taken into account in Version 2 of the arXiv listing of Ref. [1]. 

Since this is a written version of a talk given at DPF 2011, I avoid theoretical complexities and limit the 
\^ | number of equations; these can be accessed in Ref. [1]. However, I will use this report to expand on the limited 
. number of references shown in my DPF talk by including substantially more here. 

Section 2 contains a summary of the current experimental situation regarding B s mixing and B s — > f.i + fJ,~ as 
, well as SM predictions for Sr^-^+p- and AMb,- Most of the underlying theoretical formalism and derived 
relationships in this area are already well known. What is more important in comparing theory with current 
experiment is instead numerics, especially regarding the uncertainties encountered in SM predictions. 

In Section 3, we describe how the comparison of AM^ SM ' with AA/^ Expt - ) yields phenomenological bounds 
on the NP i? s -mixing contribution AM B NF \ This can, in turn, be used to constrain the magnitude of NP 
contributions to the B s — » decay mode. We implement this approach by exploring five possible NP 

models: an extra gauge boson Z' , family symmetry, R-parity violating supersymmetry, a fourth sequential 
quark generation and flavor-changing Higgs models. 
Our Summary appears in Section 4. 



2. Standard Model Analysis and Experimental Review 

The numerical predictions we cite for Sr^ 1 ^ + , AM B ^ M \ etc are based in part upon inputs displayed in 
Table I. These entries were the ones available when the analysis in Ref. [1] was being performed. Of course, 
improvements in such inputs over time will modify our results. 



Table I: Some Input Parameters 



a s (M z ) = 0.1184 ± 0.0007 [6] 


\Vu\= 0.04031™ [7] 


AM Bs = (117.0 ±0.8) x 10" 13 GeV [7] 


Ar B ./TB.=0.092±g : ™ji [7] 


mt{mt) = (163.4 ± 1.2) GeV [8] 


Jb b \J Bb s = 275 ± 13 MeV [9] 


B Bs = 1.33 ± 0.06 [9] 


f Bs = 0.2388 ± 0.0095 GeV [9] 



In Table I, the values for the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements |Vt s | and |Vtb| are taken 
from the global fit of Eq. (11.27) in Ref. [7]. Another route to \Vt s \ results in the value given in Eq. (11.13) of 
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Ref. [7], | Vi S | = 0.0387± 0.0021. In this instance, we have adopted the one in Table I because it has the smaller 
listed uncertainty. Finally, corresponding to the running mass m t (m t ) in Table I is the top quark pole mass 
tt4 po1c) = 173.1 ±1.3. 

2.1. About B s fi + n~ 

Consider first the SM prediction for Br Bs ^^+^- . A direct calculation of this quantity in the SM gives 



9 

1-4-^ 

Ml 



1/2 



where r\y is a QCD factor, Y(x t ) is an Inami-Lin function [3] and x% = m\(mt) j M"^ where mt is the running 
top-quark mass in MS renormalization. As pointed out in Ref. [4], inserting the SM expression for AMb s 
removes the |V,.*Vtb| 2 factor, thereby reducing the uncertainty in the SM branching ratio expression. It is this 
modified expression we use to obtain 1 



Br (SM) J_ AMg xpt) {G F M w m^ Y Yf 



2 

mf. 
Ml 



1/2 

= (3.3 ±0.2) x 10~ 9 , (1) 



where the largest source of uncertainty arises from P>b s followed by the dependence in So(xt) (c/. discussion 
preceding Eq. (4)) on the i-quark mass. 

In response to a request from a Session Organizer, my talk included a summary of the experimental situation 
for Br Bs ^^+ M - . It has begun to change rather substantially. Below are values taken respectively from the 
Particle" Data Group (PDG) [7], the CDF collaboration [10] and a compilation (LHC) of LHCb and CMS 
data [2], all in units of 10 -9 : 

< 47 x 10- 9 (CL = 90%) ; Brg^> + „_ = (iSt 11 ) x lO" 9 

Br £-2 +M - < 9 x 10-9 ( CL = 90% ) ; Br B™l +l >- < 11 x 10-9 ( CL - 95% ) ■ ( 2 ) 

The strongest bounds (from LHC), which have become publically available at the 2011 summer conferences, 
show just how rapidly data for the B s — > mode is being accumulated. The 'NP Window', currently 

is expected to close sometime in 2012 if B s — > n + is not observed by then [11]. As such, we are either getting 
near to encountering NP in this mode or testing the SM prediction. Either outcome is eagerly awaited. 

2.2. About AM Bs 

The experimental value [7] for AMb s , 

AMj; &pt) = (117.0 ±0.8) x 10~ 13 GeV , (3) 

is a very accurate one - the uncertainty amounts to about 0.7%. The next-to-leading (NLO) SM formula 
is arrived at from an operator product expansion of the mixing hamiltonian [12, 13]. The short-distance 
dependence in the Wilson coefficient appears in the scale- insensitive combination r)B„So(xt), where the factor 
So(xt) is another Inami-Lin function [3]. Our determination yields So(xt) — 2.319 ± 0.028. Also, we obtain 
t]b s = 0.5525 ± 0.0007 for the NLO QCD factor. From the numerical inputs discussed thus far, we find 

AM (SM) = (G F MwJ s V- tb |) 2 MgJ ^^^ 5o( _ t) = (125 ,2+13.8) x 10 -i3 GeV (4) 
which is in accord with the experimental value of Eq. (3). 



1 For another recent evaluation, see Ref. [5]. 



Proceedings of the DPF-2011 Conference, Providence, RI, August 8-13, 2011 



3 



3. New Physics Analysis 

Let us first obtain a numerical (la) bound on the New Physics contribution to AMb b - 2 We then use this to 
constrain couplings in a variety of NP models and thereby learn something about the B s — > transition. 
Accounting for NP as an additive contribution, 

AM ( l xpt) = AM { ™ } + AM^ P) , 

we have from Eqs. (3), (4), 

AM ( l P) = (-8.2±H;f) x 1(T 13 GeV . 

The error in AMj oxpt) has been included, but it is so small compared to the theoretical error in AM S (SM) as to 
be negligible. The la range for the NP contribution is thus 

AM { g s F) = (-20.9 -> +5.6) x 1(T 13 GeV . 

To proceed further without ambiguity, we would need to know the relative phase between the SM and NP 
components. Lacking this, we employ the absolute value of the largest possible number, 

|AM^ P) | < 20.9 x 10~ 13 GeV , (5) 

to constrain the NP parameters. 

Next is the issue of which model of NP to adopt. There are, in fact, a number of ways that NP can impact 
the SM: 

• Extra gauge bosons (LR models, etc) 

• Extra scalars (Multi-Higgs models, etc) 

• Extra fermions (Little Higgs models, etc) 

• Extra dimensions (Universal extra dimensions, etc) 

• Extra global symmetries (Supersymmetry, etc) 

For an analysis of D a mixing which includes as many as twenty-one NP models, see Ref. [14]. Here, however, 
we shall content ourselves with a smaller NP menu of five items. For each, we shall include a few, somewhat 
informal, introductory remarks. 

Many NP models have multidimensional parameter spaces whose complexity hinders the ability to utilize the 
constraint of Eq. (5). Several strategies come to mind for addressing this situation. Our approach is possibly 
the simplest one - to employ whatever set of reasonable assumptions (and/or physical intuitions) which allow 
us to find paths in parameter space which relate |Am]j NP ^| to Br Bs _+^+ . However, we plan to revisit this 
issue in future work. 

3.1. Single Extra Vector Boson Z' 

We should clarify what is meant by the title of this subsection. There are many Z' models, literally hundreds, 
which contain one or more of these hypothetical gauge bosons. What we have in mind is the subset of such 
models in which one Z' has a much lower mass than other NP degrees of freedom. This explains the 'single' 
extra vector boson Z 1 . In this case, the NP contribution to B s mixing arises from the Z 1 pole diagram (cf 
Fig. 1). Moreover we assume that the Z' has SM couplings to lepton pairs, thus leaving us with two unknowns, 
the Z' mass Mz> and the nondiagonal flavor coupling g Z i s j y 
The scaling with Z' parameters then goes as 

AMf oc % and Brfl + oc . ^2 



2 The possibility of utilizing Ars s is not considered here, but is addressed in Ref. [1]. 
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s " * ^ b 

Figure 1: B s Mixing via Z' exchange 



leading to 



1 TeV 



M z > 

This value will lie below the corresponding SM prediction (Sr^^^ ^_ — 3.3 x 1CT 9 ) even for a Z' mass as 
light as Mz' > 0.3 TeV. So, this class of models will not dominate the SM result. 

3.2. Family ('Horizontal') Symmetries 

The motivation for Family Symmetries (FS) is to obtain an understanding (still needed!) of fermion masses 
and fermion mixing matrices. This subject has a relatively long history, first having been actively explored in 
the late 1970's. However, it remains an area of theoretical interest up to recent times (e.g. see Ref. [15] where 
FCNC and CPV patterns are studied in the context of Family Symmetries). 

Here is a brief overview, with details left for Ref. [1]. The gauge sector in the Standard Model has a large 
global symmetry which is broken by the Higgs interaction [16]. By enlarging the Higgs sector, some subgroup 
of this symmetry can be imposed on the full SM lagrangian and the symmetry can be broken spontaneously. 
This family symmetry can be global [17] as well as gauged [18]. If the new gauge couplings are very weak or 
the gauge boson masses are large, the difference between a gauged or global symmetry is rather difficult to 
distinguish in practice [19]. In general there would be FCNC effects from both the gauge and scalar sectors. 
GHPPY analyze the gauge contribution. Consider the family gauge symmetry group SU(3)c acting on the 
three left-handed families. Spontaneous symmetry breaking will render all the gauge bosons massive. If the 
SU(3) is broken first to SU(2) before being completely broken, we may have an effective 'low' energy symmetry 
SU(2)q. This means that the gauge bosons G = {Gi} (i = 1, . . . , 3) have masses m, which are much lighter 
than those of the {G k } (k = 4, . . . , 8). 

There is a history of applications of family symmetry in which the number of unknowns becomes reduced 
to manageable proportions. The whole story is rather involved, but the following gives an impression of some 
steps. As regards the gauge boson masses {rrij}, in a simple scheme of symmetry breaking [20] one obtains 
mi = m 3 . There are in the family symmetry also four mixing matrices Ud,U u ,Ue, U v , unknown except for the 
constraints 

UlU d = V C km and UlUi = Vmnsp , 

where Vckm and Vmnsp are respectively the well-known Cabibbo-Kobayashi-Maskawa and Maki-Nakagawa- 
Sakata-Pontcorvo mixing matrices for quarks and leptons. Through reasoning given in Refs. [21], [22], it is 
somewhat possible to tame the zoo of unknown mixing parameters. 
The upshot of all this is 

AMf 1 cx and i3r ( R FS ' . oc ^ , 



where / sets the scale of the interactions between the gauge-bosons and fermions. The above relations then 
yield the result 

£r ( D FS L + .._ <0.9x 10~ 12 , 



which is tiny compared to the SM prediction. 
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Figure 2: B s — > £ + l via sneutrino exchange 



3.3. R-parity Violating Supersymmetry 

One of the models of New Physics that has a rich flavor phenomenology is R-parity violating (RPV) SUSY. 
R-parity distinguishes between particle and sparticle as 

Rp = r\3(B-L)+2S = [ +1 (particle) 
1—1 (sparticle) 

If R-parity is conserved and the initial state consists solely of ordinary matter, then intermediate and final states 
can contain only even numbers of sparticles. In this subsection, we lift this restriction by allowing the presence 
in the superpotential of terms 

Wft = \\ijkLiLjEi + KjkLiQjDt + ±% jk U?DfD° k 

where i = 1,2,3 are generation labels. We impose baryon number symmetry by setting A" to zero. Also, we 
will assume CP-conservation, so all couplings Ayft and A^- fc become real- valued. The set of {Ayfc} occur for 
couplings to leptons whereas only the {A^ fc } occur in B s mixing. For example, the lagrangian describing RPV 
SUSY contributions to B s mixing can be written as 

£fl = -Ks,2^tJ>RSL ~ A- 23 ^ H Sfl&L + h.C , 

The crucial difference between studies of RPV SUSY contributions to phenomenology of the up-quark and 
down-type quark sectors is the possibility of tree-level diagrams contributing to B s -mixing 3 and B s — > £ + £~ 
decays [23-26]. Then a process like the one in Fig. 2 dominates. Note how the initial and final states consist of 
ordinary matter whereas the intermediate state has a single sneutrino - this is clearly B s mixing due to RPV 
dynamics. 

In RPV-SUSY, the underlying transition for B s — > is s + b — > fi + + fi~ via tree-level u-squark or 

sneutrino exchange. In order to relate the rare decay to the mass difference contribution from RPV SUSY 
AMg , we need to assume that the up-squark contribution is negligible. This can be achieved in models where 
sneutrinos are much lighter than the up-type squarks, which are phenomenologically viable. Then it can be 
shown that the dependence of Br^ on the RPV parameters becomes 

Kr m or ( XmX i3A 

ar B s ^+n- « y M 2 j ■ 

Upon inserting the information from AAfg and assuming Xk23 = ^k32 , it is possible to plot the dependence of 
Brg^^+n- on Xk22 for different values of Mp i , which we present in Fig. 3. 

3.4. Fourth Quark Generation 

One of the simplest extensions of the Standard Model involves addition of the sequential fourth generation 
(often denoted as SM4) of chiral quarks, t' and b' [27-29]. The addition of a sequential fourth generation of 
quarks leads to a 4x4 CKM quark mixing matrix [30]. This implies that the parameterization of this matrix 



3 We assume that there is no strong hierarchy between the RPV SUSY couplings that favors possible box diagrams. 
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Figure 3: Branching ratio of £>r B o_ >M + M - as a function of RPV leptonic coupling \k22 and sneutrino mass Mz i = 100 GeV, 
150 GeV, and 200 GeV (solid, dashed, and dash-dotted lines). 



requires six real parameters and three phases. Besides providing new sources of CP-violation, the two additional 
phases can affect Br Bs ^^+^- due to interference effects [31]. 

There are several existing constraints on the parameters related to the fourth generation of quarks, including 
direct searches, CKM unitarity tests, and fitting precision electroweak data (S and T parameters) [32-35]. The 
latter strongly constrains the masses of the new quarks. Finally, for the sake of completeness we take note of 
a data input which became available subsequent to our DPF 2011 talk - that null results in Higgs searches by 
the LHC detectors now place the scenario of a sequential fourth quark generation 'in deep trouble' [36]. 

The relationship between AMg 8 and Br^^ + _ nas been previously studied in detail in Ref. [37]. We shall 
update their result. In SM4, the branching ratio for B s — > can be related to the experimentally-measured 4 





m t , GeV A[, s 



Figure 4: Left: branching ratio of Sr s o_ +M + M - as a function of the top-prime mass m t i for different values of the phase 
4>t' s = 0,7r/2,7r (solid, dashed, dash-dotted lines) and \l s = |Vt' s V^i,| ~ 10 -4 [39] (see also [40]). Right: branching ratio 
of Sr s o_> M + M - as a function of the CKM parameter combination X\ s with tf> t i 3 — and different values of m t i = 400 GeV 
(solid), 3 500 GeV (dashed), and 600 GeV (dash-dotted). 



4 We shall use AMj Xpt ' , as the separation of NP and SM contributions used in the rest of this paper, viz AMg B = AM^ M ' 
mL SM4 ', is not possible here due to loops with both t' and t, c, or u quarks. 
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AA4 Expt) as [37] 



Br^ M4 K _ oc AA4 Expt) ^' 2 



B s ^+H- ^ ^ IV1 B 3 | A , 

where the parameter A' is a i? s -mixing loop parameter [37], 

A' = n t S {x t ) +rjt>B%, t S (x t ') + 2rit>Rt>tSo(xt,x t >) , 

R t ' t = Vt' s V t 1 b /Vt s V t * b and the definition of the function So {x t , x t > ) can be found in Ref. [37]. The Wilson 
coefficient C{q 1 is defined as 

cli t ( P ) = cM + RftCiM 

with C 10 obtained by substituting m t > into the SM expression for Cio [38]. Our results can be found in Fig. 4. 
Referring to the recent LHC bound of Eq. (2) on Br B ^^ L +^- , the exclusion regions on the B s — > branching 
ratio in Fig. 4 lies above 1. x 1CP 8 . As one can see, the resulting branching ratios are for the most part lower 
than the LHC experimental bound. However, with values of the CKM4 matrix element X\ s = |Vt' s V t ? b | of about 

0.01, disfavored by [39], but still favored by [40], Br^ M ^\ _ can still exceed the current bound. 



3.5. Flavor-changing Neutral Current Higgs Bosons 

Many extensions of the Standard Model contain multiple scalar doublets, which increase the possibility of 
FCNC mediated by flavor non-diagonal interactions of neutral components. While many ideas exist on how 
to suppress those interactions (see, e.g. Refs. [41-43]), the ultimate test of those ideas would involve direct 
observation of scalar-mediated FCNC. The interest in multi- Higgs structures has remained fairly constant over 
the years and continues unabated to this day (e.g. see Ref. [5]). 

A generic interaction hamiltonian of this type is 



H 



H 



V2 



'V2 



>^23 b R s L 



+ i\§2 t>LS 



R + I 



H.c. 



where <f>° and a° represent the lightest scalar and pseudoscalar states respectively and ellipses stand for the 
terms containing still heavier states whose contributions to AMb s and Br Bs ^.^+^- will be suppressed. We take 
all couplings and Xfj to be real- valued. The superscripts D and E on these refer respectively to couplings of 
d-type quarks and of charged leptons. To proceed, we need to distinguish two cases: the lightest FCNC Higgs 
particle as a scalar (0°) or pseudoscalar (a ). 
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Figure 5: Branching ratio of Br B o_ >tt + tt - as a function of pseudoscalar Higgs mass Mo. Left: A22 = 1,0.5,0.1 (solid, 
dashed, dash-dotted lines). Right: \f 2 = 0.1,0.05,0.01 (solid, dashed, dash-dotted lines). 
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3.5.1 . Light scalar FCNC Higgs 

The case of relatively light scalar Higgs state is quite common, arising most often in Type-Ill two-Higgs doublet 
models (models without natural flavor conservation) [44-46] . Although the FCNC Higgs model does contribute 
to AMb,, it does not contribute to B s — > at tree level. Any nonzero contribution to B s — > decay 

must be produced at one- loop level [47]. 

3.5.2. Light pseudoscalar FCNC Higgs 

The case of a lightest pseudoscalar Higgs state can occur in the non-minimal supersymmetric standard model 
(NMSSM) [48-51] or related models [52]. In NMSSM, a singlet pseudoscalar is introduced to dynamically solve 
the /i problem. The resulting pseudoscalar can have a mass as light as tens of GeV. This does not mean, 
however, that it necessarily gives the dominant contribution to both B s mixing and the B s — > decay 
rate since there can be loop contributions from other Higgs states. However, here we work in the region of the 
parameter space where it does. Calculation reveals 

AM (a) oc (*™\ and Sr (a) oc ( X * 2 

In Br^g t+l _ , the dependence on \® 2 can be eliminated by using AMg' . The unknown factors then enter in 
the combination A^/M a . In Fig. 5, we plot the dependence on M a for different values of Af 2 - 

4. Summary 

Our talk consisted of two main parts, the (updated) SM evaluations of Ref. [1] and the issue of NP contribu- 
tions to B s — > Wc discuss each of these in turn. 



4.1. Update of SM Evaluations 

As regards B s mixing, we have 

AM (Expt) = (H7.0 ±0.8) x 10~ 13 GcV and AM ( B SM) = (125.2±H;|) x 1(T 13 GeV . 

The uncertainty in the SM result is seen to be about 16 times larger than that in the experimental listing. It 
arises mainly from the factors f B Bs s and |Vt s | 2 in Eq. (4) and is roughly equally shared between them. This 

large theoretical uncertainy in AAff"' hinders the study of additive NP contributions. 

As regards the B s — > branching ratio, upon using as input (c/Eq. (1)), we obtain 

Br^ + _ = (3.3 ± 0.2) x 10~ 9 . 

We find the main uncertainty to be from Bb s and roughly half as much from the implicit f-quark mass depen- 
dence in So(xt). 

4.2. NP Contributions to B s — y 

The GHPPY approach to B s -4- is to use 

AMg xpt) =AMg M) +AM^ P) 

to bound |AM^ P '| and thus to constrain NP parameters. In this talk, I described the results of studying five 

NP models, with a sixth in preparaion. For two of them (the cases of a single Z' and of family symmetry), we 
conclude that the contribution to B s — > n + lies below the SM prediction. For the others, Figs. 3-5 depict 
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the regions of exclusion in their respective parameter spaces. We repeat our earlier comment that the exclusion 
regions in earlier versions [1] of our Figures are based on the PDG listing 

Br £^2+ M - < 47 x 10-9 (CL = 90%) 

and have been updated here to reflect the combined LHCb and CMS bounds, 

£ r (LHC) f 9 x 1(T 9 (CL = 90%) 

r B 3 ^+»- j n x 10 _ 9 ( CL = 95% j 

As a final comment, we point out that the size of the SM contribution to D° — > [i + [i~ [53] is quite feeble 
compared to that in B s — > \i + ■ It would appear that the D° — > decay mode could well be a fruitful 

arena to search for New Physics! 
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